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Abstract: Based on the probability generating function of stuttering Poisson distribution (SPD), this 
paper considers some equivalent propositions of SPD. From this, we show that some distributions in 
the application of non-life insurance actuarial science are SPD, such as negative binomial distribution, 
compound Poisson distribution etc.. By weakening condition of equivalent propositions of SPD, we 
define the generalized SPD and prove that any non-negative discrete random variables X with P{X = 
0} > 0.5 obey generalized SPD. We consider cumulant estimation of generalized SPD's parameters. 
As an application, we use SPD with four parameters (4-th SPD) to fit auto insurance claim data. The 
fitting results show that 4-th SPD is more accurate than negative binomial and Poisson distribution. 
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1 Introduction 

Stuttering Poisson distribution (simply write SPD) is a non-negative discrete compound Poisson dis- 
tribution (see[16],[12]), which has the feature that two or more events occur in a very short time (arrive 
in group or batches). For example, a man may claim for double or much compensation because he 
has two or more insurance policy from the insurance company. In inventory management, a customer 
may buy more than one goods of same kind. 

Definition 1.1 For a stochastic process {£(i),i > 0}, let 

P n (t) = P{£(t)=n|£(0)=0}. 

Similarly to some properties of Poisson process [14], stuttering Poisson process £(t) satisfies the fol- 
lowing properties: 

Property 1: £(0) = 0. 

Property 2: £(t) has independent increments (i.e., the numbers of events that occur in disjoint 
time intervals are independent) and stationary increments ( i.e., the distribution of the number of 
events that occur in any interval of time depends only on the length of the time interval). 

Property 3: Pi(At) = Ac^At + o(Ai), (0 < a* < 1, i = 1, 2, ■ ■ ■ ,r), where 1 < r < +oo. The r 
equal to finite or infinite in the following part of this paper. 
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Property 3 implies that stuttering Poisson process satisfies that the probability of two or more 
events occur at a very short time interval is non-zero. The probability is directly proportional to the 
length of time interval. When i = 0, according to the properties of independent increments , stationary 
increments and Chapman-Kolmogoroff equations, we have 

P {t + At) = P (t)P {At) 

implying that 

P {At)=e~ XAt =l - XAt + o(At), (A > 0) 
Add terms from Po(At) to P r (At), that is 

r r 

1 = 5^ Pi (At) = 1 - XAt + o(At) + AaiAt + o(At), 

i=l i=l 

hence we obtain Ya=i a i = 1- 

Definition 1.2 If satisfies Property 1 ,2 and 3, then we say that £(t) obeys the r-th stuttering 
Poisson distribution. Denote 

~ SP(ai\t, ■ ■ ■ ,a r Xt) 

with parameters (a±Xt, ■ ■ ■ , a r Xt) € W , (a r 7^ 0). 

Similar to Poisson law of small numbers, SPD can be deduced from the limiting distribution 
of multinomial distribution [19]. Let pi = ^r, then the probability generating fubction (PGF) of 
stuttering Poisson distribution 

P(s) = lim [(1 - J>) + = lim [1 + -E°*f = e A *£ ai(-, " 1) . 

Npi=cti\t * — ' * — ' N^-oo iv * — ' z — ' 

jy^OQ i=l i=l i=l i=l 

where on(i = 1, 2, ■ ■ ■ ) is probability density of an positive discrete distribution. With P n {t) = P< J/ ^ 
and Fa di Bruno fomula [15] 



^ — 1 J.. _l Li- J Z. * Z. ^ 



= 1 fciH hfe u +---fc„=i,fc u eN 

1-feiH hwfc-tiH hrafcn=n 
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we have 

P„(t) = KAt + ...+ ^ aUs! • • • ( } + ' " + ^!^ ]e ' (A > 0) - (1) 

fclH hfc u + ---fc„=i,fc ll GN 

1-felH huk u -\ \-nk n =n 

Other methods to prove the expression of stuttering Poisson distribution can be obtained by the 
limiting distribution of multinomial distribution [19], system of differential equations [17], system of 
functional equation [9]. The name compound Poisson was used by W. Feller [15] and R. M. Adelson 

At E ^(s'-l) 

[10] to discuss distribution which PGF is e i=1 . When r = 1, SPD degenerates to Poisson 

distribution. When r > 2, we call it non-degenerative SPD. When r = 2, CD. Kemp and AW. Kemp 
[2] named it Hermite distribution owing to the PGF can be expanded in terms of Hermite polynomial. 
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When r = 3(4), Y. C. Patel [17] said it triple (quadruple) stuttering Poisson distributions. H. P. 
Galliher,et al. [7] considered the demands £(t) obey SPD with parameters 

((1 — a)X, (1 — a)a\, (1 — a)a 2 \, • • • ) 

of geometric distribution in inventory management theory and first name it stuttering Poisson. T.S. 
Moothathu, C.S. Kumar [12] considered SPD with parameters of binomial distribution. For more 
application in inventory management, see [3], [6], [8]. In queue theory, A.Kuczura [1] considered requests 
arrive in group or batches with constant service times. R. Mitchell [11] showed that SPD is a more 
exact model than the Poisson model to fit observed demand form actual historical data of several U.S. 
Air Force bases. D. E. Giles [4] used Hermite distribution to fit data for the number of banking and 
currency crises in IMF-member countries. 



2 Equivalent Propositions of SPD 

Let P(k,p) be a family of non-negative discrete distributions, which is closed under convolution oper- 
ation, where p (denoting the mean value of the distribution) runs over all non-negative real numbers. 
L. Janossy [9] showed that P(k,p) is SPD, The equivalent proposition of SPD is deduced from PGF 
and it help us to consider the generalized stuttering Poisson distribution. 

Theorem 2.1 For a discrete random variable X = = 0, 1, 2, • • • ) with P(s) = YH=oPi sl (\ s \ — 
1). Then, taking logarithm of PGF and expand it to a power series 

r r 

g x {s) = InC^piS 1 ) = ^bis\ \s\ < 1 
i=0 i=0 

with J2i=o h = A < oo, {pi > 0), where gx(s) is cumulants generation function of a discrete random 
variable. Then, if and only if the discrete random variable obey SPD. 

Proof: Sufficiency. With YH=o^ ist = m (Si=o^ iS *) an< ^ S[=o ^ = A < oo, we know that gx(s) 
is absolutely convergent in |s| < 1. Hence Ya=o h = m (Si=o^) = 0- Let ^* = a *A> it yield to 

P(z) = e !=1 *=i = e i =! 

Set X't = A, then X ~ SP(Xai, Xct2, ■ ■ ■ )• 

Necessity. The parameters of SPD satisfy Yli=o = A < oo. 

Example 2.1 Negative binomial distribution (NBD): P k = ^ ^ ^jp r (p~ l) k ,(p € (0,1), k = 

1,2, • • • ). The PGF of NBD is [ 1 _ (1 p _ p) J r . The logarithm of PGF is rlnp + £ r -^ L s\ (\s\ < 1), 

' i=i 

then we have 

r(l — pf 1 — p 
rmp+> <r(mpH ). 

h 1 p 

So the NBD is equivalent to SPD with parameters (rq, ■ ■ ■ , ■ ■ ■ ). This reveals the essential 
properties of NBD. NBD is an important distribution in automobile insurance, which is a mixture 
distribution of Poisson distribution and logarithmic distribution [12]. 
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W. Feller [15] name that the discrete random variable X(i.i.d., P{X = i} = bi,i = 0, 1,2, •••) 
whose sums Yn = YliLi X% with N ~ -P(Ai) is compound Poisson distribution (the accurate name 
should be discrete compound Poisson distribution or SPD). And PGF is 

r 

E(E{ s Yn \N = n))=E{[G(s)} n )=e x ^~^ =e £ ^ . 

SPD have infinitely divisible property via its PGF. A nonnegative discrete distribution is called 
infinitely divisible if for any n > 1, its PGF can be represented as the n-th power of some other PGF. 

Thus a SPD with PGF e i=1 , then this SPD can be represented as the n-th power of the other 

<*(.'-U %±°iV-i). _ mXt arM , 



n 



PGF e <=i . e i=1 is SPD with parameters (^f 

Next, we obtain "compound compound Poisson" sums is also SPD in Theorem 2.2 . 

Theorem 2.2 X(P{X = i} = bi,i = 0, 1,2, •••) are independent identically distributed i.i.d. 
random variables. When N ~ SP(a\\t, ••• , a r \t), "compound compound Poisson" sums Yn = 
J2iLi Xi is SPD. 

Proof: Gx(s) = Yll=o bi s \(\s\ < 1)> and using double conditional expectation of s Yn 

v v , A Y, ai[G(s)*-l] 

/V JV ( S ) = J B( S ^) = ^(E( S ^|iV = n))=^([Gx( S )n=e *=i 
Noticed that < Si=o = 1> hence 

r r 

A ^ ^[Gxisf - 1]| < A ^OilGxWI* + A < 2A. 
i=i i=i 

So Py Ar (s) is absolutely convergent in \s\ < 1, 

r r r r 

A - 1] = \J2 «i(6o + M + + • • • + MT - A = ^ Aqs 1 - Ao^l - 

j=0 i=l i=l j'=l 

where q are derived from multinomial expand. We need't to have the accurate expression of q. Let 
£^ =1 ci!j(l — b 3 )=c, then 

AEci( a *-l) 

Noticed that q > and ^[=i Aq = c < oo, Yn is SPD by using Theorem 2.1 . 

In Theorem 2.2 , a discrete compound Poisson distribution (or SPD) is a special case of discrete 
compound SPD (compound compound Poisson distribution) when N ~ P(Xt). Theorem 2.2 show 
that SPD is infinitely divisible distribution. It concludes that "compound ••• compound Poisson 
distribution" is SPD. In practice, many claims may be from superimposed events. That explains 
why some distributions in non-life insurance are equivalent to SPD. Besides SPD, other generalized 
Poisson model has wide applications in non-life insurance actuarial model and risk model such as 
mixed Poisson process [14], [17] and doubly stochastic Poisson processes [17]. 



3 Generalized Stuttering Poisson Distribution 

In this chapter, general stuttering Poisson is defined by weakening conditions of Theorem 2.1 . 
Furthermore, we show that any non-negative discrete random variable X with P{X = 0} > 0.5 obey 
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generalized SPD. L. Janossy [9] used independent increments, stationary increments and Chapman- 
Kolmogoroff equations to construction the system of functional equation 

i 

m +At) = Y, p k (&m_ k (t),(i = o, i, • • • ) (2) 

k=0 

Solving (2) from one to one will deduce to (1). For example, when i = 0, we have Po(t + At) = 
P (t)P (At), the solution of P (t) is P (t)=e~ xt . 

When i = 1, we have Pi(t + At) = P (t)Pi(At) + Pi(t)P (At), the solution is Pi(t) = aite~ A< . 

When i = 2, ■ ■ ■ , by the system of functional equation (2) we have 

P B (t) = [a B At + --- + £ ^JL - - - (At)- + • • • + -iLJ-] e -M (A > 0), (3) 

fciH hfc u +---/c n =i,fc u eN 

1-feiH l-tife u H hrafc n =« 

where (A > 0, aij G R, n = 0, 1, 2, • • • , r). Janossy [9] proved (3) by mathematical induction. 
There are no nonnegative restriction in (3), implies 

Pi (At) = \ ai At + o(At), (i = 1, 2, • • • ) 

we obtain a i = 1- 

Since ccj are not necessarily to be nonnegative, now we suppose that ojj may take negative value 
and satisfy Yli^i ccj = 1 to have a new distribution family. 

Definition 3.1 Generalized stuttering Poisson distribution (GSPD): For a discrete random vari- 
able X(X = i, % = 0, 1, 2, ■ ■ ■ ) , the form of PGF is e^ita*' -1 ) and satisfies A > 0, J^Zi a i = L 
When aj = 0, if «j > r + 1. We name it r-th generalized stuttering Poisson distribution (GSPD). 

It is obvious that SPD is a subfamily of GSPD. 

Theorem 3.1 Any non-negative discrete random variable X(P{X = i} = bi, i = 0, 1, 2, • • • ) with 
P{X = 0} > 0.5 obey GSPD. 

oo 

Proof: The PGF of X is P(s) = g pt i s i (|s| < 1). Let Q(s) = + ^s 2 + • • • , expand ln[P(s)] 

i=0 

with Taylor series at zero, that is 

oo , _. sj—1 CO 

ln[P(s)] = ln[l + Q(s)\ + lnp = ln Po + 1 , J Q»(s) = to*- 

i=l «=o 

Noticed that |Q(s)| < Pq'Hpi + f>2 + • • • )■ If Po 1 (pi + P2 + ■ ■ ■ ) < 1, then we have X)i=i ^ s 
absolutely convergent in |s| < 1. Solving the inequality, we get po > 0.5. On the other hand 

oo 

^2 bi = ln[P(l)] - lnpo = - hipo, 
i=i 

judging from Definition 3.1, X obey GSPD. 

The zero-inflated model is a random event containing excess zeros count data in the unit time [5],. 
For example, the number of claims to an insurance company is almost zero, otherwise substantial losses 
to make insurance company bankrupt. Zero-inflated model have too much count data that Poisson 
distribution is disable to forecast. It is obvious that condition in Theorem 3.1 satisfies zero-inflated 
model for P{X = 0} > 0.5. 
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Example 3.1 Theorem 3.1 show that the probability of zero occurrence is more than positive 
occurrences. For example, the cumulant of Bernoulli distribution P(s) = p + (1 — p)s,(p > 0.5) is 

Inb + (1 - P)s] = - hip + Iz£ fl - * (Iz£) V+ .... 

p 2 p 3 p 

Let t = 1, thus the other form of PDF is Pi(At) = =± (^y-J At + o(At), (i = 1,2, •••). 

Intuitive explanation: The possibility of two or more events at a given instant is permitted with 
negative and positive probability in turn. Have a "negative probability" in an instant can be considered 
that demands leave and arrive alternatively in group or superposition. Leaving and arriving with same 
probability may be seen as zero probability. 



4 Statistic of Generalized Stuttering Poisson Distribution 

Cumulants n n , moments m n and central moments c n of GSPD are deduced from probability generat- 
ing function and moment generating function. It also can use in SPD. A discrete random variable 
X(P{X = i} = bi,i = 0, 1,2, • • • ) with PGF P(s) = J2'iZoPi s \ (\ s \ < 1)> where moment generating 
function is Mx(s) = P(e s ). Expanding e sX with Taylor series at zero, we have 

OO j — i , rTI oo 



Mx(s) = E(E (£*>!) = E ^ i E ( |.| < i) 

n n n=0 

Definition 4.1 Cumulants generating function of a random variable is 

oo n 

g x (s) = HM x (s)) = J2^, (N < 1) 



?1=0 



where coefficients n n (n = 0, 1, 2, • • • ) is n-th cumulants. It is explicit that ko = 0. 
Theorem 4.1 ~ GSP(a>iXt, ■ ■ ■ ,a r Xt), then the n-th cumulant of is 

r 

n n = ^2ai\ti n . (4) 



i=0 



Proof Cumulants generating function of GSPD is gx(s) = Yli=o ctiXt(e ls — 1), Expanding e with 
Taylor series at zero, that is 

r OO j OO y' OO ,, y- 

E oa\t{e is - 1) = -Xt + At[ ai (l + E £ [) + «2(1 + E i ^ L ) + • • • + + E 

i=l j=l J ' j=l J ' j=l ■ 7 ' 

r r ^ r ^ 

= (E <*iAiz)s + (E «i^ 2 )— + ■ ■ ■ + (E «^)- + • • • . 

i=l i=l i=l 

Comparing coefficients of t n , we have n n = Y^i=o otiXti n . When r — > oo, k„ may divergence. 

Theorem 4.2 If £(t) ~ GSP(aiXt, - - - , a r At), then the recursion formula of nth moments m n is 

= E ( r ) K n+i-jmj = E ( " ) E «A^ n+1 "'')^, (ki = mi = E OiAti). (5) 

j=0 ^ ' j=0 ^ ' i=l i=l 
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Proof Expanding ln[Mx(s)] with Taylor series at zero, that is 

hm x (s)] = hi + e = E ! Jr n - 2 1 (E > n ) 2 + ■ ■ ■ + ^-(E ^ + ■ ■ ■ 

n=l n=l n=l n=l 

777,2 — m? o 777,3 — 377717772 + 2777? , 777-4 — 477737771 — 3?77o + 12777 2 777? — 6mf a 

= m ^ s+ 2! 5 + 3! " + 4! ' + 

By the definition of cumulant, we have 

r r 
«1 = 777i = EX = OiiXU => 777i = a i^i, 

i=l 1=1 

r r r 

k 2 = 777 2 - ml = E(X - EX) 2 = E aiXti 2 => m 2 = E ««-^ 2 + (E a ^) 2 ' 

i=i j=i j=i 

r 

^3 = — 377717772 + 2?nf = — EX) 3 = CtiXU 3 =4> 7773 = K3 + 377717772 — 2jT7,f , 

1=1 

K4 = m 4 - 477717773 - 3m| + 12777^7772 ~ 6777^ = E(X - EX) 4 - 3[E(X - EX) 2 ] 2 , ■ ■ ■ 

By taking the derivative of both side of Mx(s)=e 9x ^ again and again with respect to s, using Leibniz 
formula we obtain 



i=0 ^ ' 



Substitute s = to the n-th derivatives of gx{s) and Mx(s), hence M^\s) = m n and g^' '(t) 
K n . Substituting cumulants and moments into (6), we have 



s=0 



777 n+ i = E ( " ) Kn+l-jTUj = E ( " ) E ^ ^""^ ) m 3 ■ ( 7 ) 



The relationship between cumulants and moment is 

n— 1 

/ 77 
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= "7n - E ( " 1 ) K n-jmj, («i = 777i) (8) 

= f(m 1 ,m 2 , ■ ■ ■ ,m n ). (9) 

Remark Higher cumulants (77 > 4) are different to central moment when 77 is more than 4. 
Ar guing from K n — Jpj- ln[iWx(s)] and Fa di Bruno formula, it can deduce to K n too. 

Example 4.2 An alternative approach to compute the mean and variance of Bernoulli distribution 



Kn = y ai M n = y ^—^-(- in P )i n = y 



1=1 



^— ' —ip l \np ^ p 

1=1 i=i 



00 1 i 00 1 i 

EX=Y,(—) =1-P,DX = £<(—) =P"P 2 
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Theorem 4.3 If £(t) ~ GSP(ai\t, • • • , a r Xt), then the recursion formula of nth moments c n is 

Cn+l = ( " ) K n+l-jCj' ( C = 1, Cl = 1, < = { ^ ^ 

Proof Expanding ln[Mx-£x(s)] with Taylor series at zero and comparing the coefficients of 
cumulant generating function, we obtain 

r r r r 

= 0, K* 2 = ^2 ctiXti 2 , k% = ®iXti 3 , k\ = ctiXti 4 + 3(^^ aiXti 2 ) 2 , ■ ■ ■ 



i=i 



i=l 



i=l 



i=l 



Since 



i=0 



substitute s = to the n-th derivatives of gx-Ex(s) and Mx-gx(s), hence M^ £X (s 



( n ) / \ 

c n, 9x-EX\ S ) 

(10). 



s=0 



s=0 



k* . Similarly to the proof in Theorem 4.2 , displacing K n with «*, we have 



5 Cumulant Estimation of Generalized Stuttering Poisson Distri- 
bution 

Y. C. Patel [18] gave moment estimator of the parameters of Hermite distribution. Y.C.Patel [19] 
estimates the parameters of the trip and quadruple stuttering Poisson distributions with maximum 
likelihood estimation moment estimation, and mixed moment estimation of the parameter. Use Same- 
pie moments rh\ and central moments C2 and £3, we have 





2 
2 




(11) 



where Ecu = on + 0(±) , (i = 1, 2, 3). 

When n > 4, From the computing in Theorem 4.2, Higher cumulants are different to central 
moment or moment. Central moment or moment are nonlinear combination of oti{i = 1, 2, • • • ). Thus 
it is difficult to estimate the parameters by using central moment or moment estimation. Theorem 
4.2 implies that K n is linear combination of cti(i = 1, 2, ■ ■ ■ ). 

Therefore, firstly, we use sample moment K n to calculate in (4). From (4), when r < 00, K n is 
convergence. Secondly, by solving the following system of linear equations of = 1,2, ■■■ , n) by 
means of 0-th to n-th cumulants formula 



/ £ \ 



/ 1 





1 \ / -1 \ 



n 



\k n ) \ 1" 



(12) 
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where the coefficients matrix in (12) is invertible Vandermonde matrix. By solving the linear system 
equation we have 

n 

&i = ^2bijkj,(i = 1,2, ••• ,n) (13) 
Assuming the samples to the power of n, Q(i = 1, 2, • • • , m) are i.i.d.. Let samples n-th moment 

m 

A n = j arguing from Khintchine's law of large numbers, we obtain lim P{\A n — m n \ < e} = 1 

n— >oo 

for all e > 0. From relationship between cumulants and moment (4), for all e > 0. f(xi, x 2 , ■ ■ ■ , x n ) is 
continuous of several variables, it implies 

lim P{\K n -k n \ <e}= lim P{\f(mi, m 2 , ■ ■ ■ ,m n ) - f(A 1 ,A 2 , - ■ ■ ,A n )\<e} = l (14) 
Thus we prove cumulant estimation is consistent estimate by the linear relation in (13). 

6 Applications 

R. M. Adelson [17] put forward the recursion formula of SPD's probability density function by using 
Leibniz formula 

P j+1 (t) = — L_[ ai AP,(t) + 2a 2 XPj-i(t) + • • • + (j + l)a j+1 XP (t)],P (t) = e - xt , (15) 

(15) avoid tediously computing the sum of much index in (3) by recursion relation. There are no 
nonnegative restriction to Oj(i = 1, 2, • • • ), so (15) can be used in GSPD. 

Now we use SPD fitting according to auto insurance claims data (the car insurance claims data of 
the following table 1 from [17]), and then we compare the goodness of fit with some other distributions 

From data in Table 1, Total insurance policies are n = 106974. The probability of zero claim 
policies is far greater than 0.5. Obviously it is zero-inflated data. The number of the insurance policy 
of i-th is Xj(j = 1, 2, • • • , 106974), so the mean value and the 2-th and the 3-th central moment of the 
insurance policy claims rate is 

mi = 0.1010806364, c 2 = 0.1074468102, c 3 = 0.1216468798. 

According to the (8) and (11), we have 

ai = 0.97255, a 2 = 0.02496, d 3 = 0.00249. 

Thus we can infer that the probability of claims of customer who buy two copies of the same insurance 
is only 2.496%, and three copies of the same insurance is only 0.249%. Employing recursion relation 
(15) and (8), we obtain p{. And then figure out npi. Analogously, consider quadruple SPD fitting, in 
this case we have 

ai = 0.97151, a 2 = 0.02703, a 3 = 0.00112,a 4 = 0.00034. 

In Table 1, we assume the data come from Poisson distribution, triple SPD (by recursion formula 
(15)), quadruple SPD (by recursion formula (15)), negative binomial distribution, respectively, and 
then estimate the probability of different numbers of claims. 
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Table 1: The comparison of auto insurance claims data of different distributions fitting effect (moment 
estimation or cumulant estimation) 



i 


Methods 





1 


2 


3 


4 


i > 4 


Vi 


Observed frequency 


96978 


9240 


704 


43 
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npi 


Estimate by Poisson 


96689.5 


9773.5 


494.5 


16.6 


0.4 





npi 


Estimate by triple SPD 


96974.1 


9256.0 


679.2 


60.4 


4.0 


0.3 


npi 


Estimate by quadruple SPD 


96977.3 


9243.2 


697.6 


49.1 


6.1 


0.7 


npi 


Estimate by NPD 


96985.4 


9222.5 


711.7 


50.7 


3.5 


0.2 



4 2 

Constructing test statistical:^ = J2 — n, by calculating, we get 

i=0 

r) PD = 345.1250, m-SPD = 12.5786, m-SPD = 2.8963, tjnbd = 10.1294. 

From Pearson's chi-squared test theory, in one hand, xl(O.Ol) = 12.277, given significant level of 
0.1 we accept that claims data obey quadruple SPD or negative binomial distribution. In the other 
hand, x|(0.5) = 3.357, given significant level of 0.5 we accept that claims data obey quadruple SPD 
fitting. It is thus clear that quadruple SPD fitting effect is better than that of negative binomial 
distribution fitting effect, and NBD is better than triple SPD. The goodness of Poisson distribution 
model is worst in those four distributions. 
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